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The ground state of the two-leg Heisenberg ladder is identified as an RVB type spin liquid, which 
is generated by Gutzwiller projection of tight-binding bands with flux 7r per plaquet. Explicit trial 
wave functions for the magnon and hole excitations are formulated in terms of spinons and holons. 



Of 
< 



i 

CO 



i 

C 

O 
o 



> 

(N 
O 

o 
o 

OO 
c3 



c 

o 
o 



X 
S3 



PACS numbers: 71.10.-w, 75.10.-b, 75.90,+w, 74.72.-h 



1. Itinerant antifcrromagncts confined to coupled 
chains, or t-J ladders, have enjoyed enormous popu- 
larity over the past few years They provide the 
simplest example of a generic spin liquid in dimensions 
greater than one, and the only example thereof which 
is presently fully amenable to numerical methods. (It 
is furthermore widely believed that they constitute the 
first step towards understanding the two-dimensional t-J 
model starting from one dimension, but I rather believe 
the two-leg ladder to be just a special case.) These mod- 
els are approximately realized in (VO)2P207, SrCu203, 
and Sri4_ a; Ca 2 ;Cu2404i, and hence accessible to experi- 
ment. 

The t-J Hamiltonian for the ladder is given by 
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where (tij, Jy) = (t, J) if i and j are nearest neighbors 
along one of the chains, and (tj_, J±) if they are nearest 
neighbors across the rungs; each pair (ij) is summed over 
twice and no doubly occupied sites are allowed. 

One of the most striking features of the two-leg t-J 
ladder is the persistence of a spin gap A » J±/2 in the 
weak coupling limit J± <C J. (For sufficiently strong 
couplings Ji > J, the system can be described by a 
perturbative expansion around the strong coupling limit 
consisting of singlets across the rungs Jl(]] , which yields a 
spin gap A « J± — J + i J 2 / J±; a weak coupling expan- 
sion starting from decoupled chains, however, is not pos- 
sible, as the individual spin chains are quantum critical 
in the sense that the tiniest perturbation can change the 
universality class.) In this letter, I will formulate a mi- 
croscopic theory of the two-leg t—J ladder, which is uni- 
versally valid at all ratios Jj_/J, in terms of explicit spin 
liquid trial wave functions for the ground state, magnon 
(spinon-spinon bound state) and the hole (holon-spinon 
bound state) excitations. 

2. The trial wave function for the ground state of the 
Heisenberg ladder, the t-J ladder without any holes, is 
constructed as follows. Consider a tight binding ladder 
with flux 7T per plaquet, and hopping terms of magnitude 
t along the chains, and t±_ across the rungs. In the gauge 
depicted in figure [l], we write the single particle Bloch 
states 
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FIG. 1. Flux band structure of a tight binding ladder with 
flux n per plaquet for t±/t = (doted lines) and tj_/t = 1 
(solid lines). The energy gap is given by 2t±. 
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where the u q (j) are strictly periodic in both real and 
momentum space and obey 
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with m = tj_/2t. Since (^) is a two-dimensional Dirac 
equation, i.e. is diagonal, we immediately obtain the 
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The coupling between the tight-binding chains hence in- 
duces an energy gap of magnitude 2t±. We now fill the 
lower band twice, once with up-spin electrons, and once 
with down-spin electrons; the resulting Slater determi- 
nant IV'sd) is obviously a spin singlet. The spin liq- 
uid trial wave function for the Heisenberg ladder with 
Jj_/J = t±/t is obtained by eliminating all the doubly 
occupied sites via Gutzwiller projection, 
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Since the Gutzwiller projector Vq commutes with the to- 
tal spin operator, Atrial) is also a singlet. This trial wave 
function is as accurate an approximation as the Haldane- 
Shastry state jjT[ for the one-dimensional Heisenberg 
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TABLE I. Energy expectation values and nearest neighbor 
spin correlations for the spin liquid trial wave functions in 
comparison with the exact ground states of a 2 x 10 Heisenberg 
ladder with periodic boundary conditions, as well as over- 
laps between trial wave functions and exact ground states. 
Throughout this article, all energies quoted are in units of 
max(Jx,J). The boundary phase before Gutzwiller projec- 
tion has been 0. 



chain in the weak coupling limit J±/J = 0, and exact 
in the strong coupling limit Jj_/J — > oo; the approxima- 
tion has its worst point at isotropic coupling (see table 

i mm- 

3. There are essentially two ways of constructing spinon 
and holon excitations for spin liquids (they are obtained 
from each other by annihilating or creating an electron on 
the spinon or holon site). The first one is Anderson's pro- 
jection technique pj : inhomogenities in both spin and 
charge created before Gutzwiller projection yield inho- 
mogenities in spin only after projection. Anderson writes 
a state with two spinons localized at sites i and j 
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In the case of the ladder, however, the spinons are not free 
particles, but bound into pairs by a linear confinement 
force Jig]. To obtain the magnon trial wave function, 
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a hence nontrivial internal wave function (f>i,j{k) for the 
spinon-spinon bound state is required. (In the one- 
dimensional spin chain, by contrast, the spinons interact 
only weakly, and the internal spinon wave function can 
be approximated by plane waves.) 

It is therefore expedient to use the second method, 
which has been successful in describing the fractionally 
charged solitons in polyacytelene Jig] : Rokhsar JlTj con- 
structs elementary excitations of spin liquids via localized 
midgap states, which are either occupied by a single elec- 
tron (spinon) or left unoccupied (holon). The topology 
of the ladder dictates that midgap states can only be 
created in pairs, which implies that we automatically ob- 
tain spinon-spinon (or holon) bound states rather then 
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FIG. 2. Magnetic tight-binding configuration for a holon of 
the chiral spin liquid (a) as proposed by Rokhsar (cutting all 
the links to a given site and adjusting the flux according to 
Rokhsar's loop rules generates a stationary holon or holetto) 
or (b) by combining this flux adjustment with Anderson's pro- 
jection technique (only the phases of the hopping parameters 
around the holon site are adjusted). 



isolated spinons (or holons). The fact that the energy 
required to create two singly occupied midgap states is 
proportional to t± suggests a spin gap proportional to J±. 
These general observations, however, leave us still with 
a large number of possible choices for the midgap states; 
most constructions yield satisfactory magnon, but only 
very few acceptable hole trial wave functions. To identify 
those, let us step back and take a broader view. 

4. The spin liquid proposed above is, in fact, a special 
case of the Kalmeyer-Laughlin chiral spin liquid jl8| , ob- 
tained by imposing a periodic boundary condition with 
a periodicity of only two lattice spacings in y-direction. 
This chiral spin liquid may be generated from a tight- 
binding lattice with flux ir per plaquet and hopping mag- 
nitudes t and ^t± in x- and y-direction, respectively; the 
P and T violating diagonal hopping elements, which are 
otherwise required to open an energy gap, cancel due to 
the boundary condition. 

Spinons and holons for the chiral spin liquid may be 
constructed via Anderson's method or via midgap states; 
Rokhsar creates a midgap state in the flux band structure 
before projection by cutting all the links to a given site 
and adjusting the flux according to his loop rules, which 
require that the kinetic energy on the loops around each 
plaquet is minimal (see figure ||a) . The resulting holon 
is not nearly as mobile as Anderson's, but optimal with 
regard to the magnetic energy; it adequately describes 
stationary charge excitations. I call it a stationary holon 
or holetto. To obtain the generic and mobile holon, wc 
create a midgap state by adjusting the flux according 
to Rokhsar's procedure (i.e. we create a defect of flux n 
around the holon site) without cutting any links (i.e. we 
adjust the hopping phases without adjusting the magni- 
tudes), and then project such that this site is unoccupied 
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FIG. 3. Magnetic tight-binding configurations and the cor- 
responding density of states for the ladder with a spinon 
bound to (a) a stationary holetto or (b) a mobile holon be- 
fore Gutzwiller projection. Only the latter flux configuration 
violates P and T. 



(see figure gb). This holon is equivalent to Anderson's in 
the case of the chiral spin liquid, but more generally ap- 
plicable. 

5. The flux configurations used to construct holetto- 
spinon and holon-spinon bound states for the ladder are 
shown in figure ^a and b. In the case of the holetto, we 
cut all the links to a given site; as the topology of the 
ladder does not provide a context for a flux adjustment 
around this site, we obtain a second midgap state, and 
hence a spinon, localized nearby. This trial wave function 
describes a stationary hole. To construct a mobile hole, 
we create the midgap states by only adjusting the flux, 
and project such that the holon site is unoccupied; as we 
are creating two rather than one midgap state, we remove 
flux 7r/2 from each neighboring plaquet fujfl . The flux 
configuration now violates P and T, and the holon-spinon 
bound state carries a chirality quantum number, which is 
+ for the configuration shown in figure ||b, and — for its 
complex conjugate; states of opposite chiralities map into 
each other under P or T. The final trial wave functions 
for the hole is a linear superposition of the holon-spinon 
bound states of both chiralities at each momentum, 
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where Tj is the holon cordinate and a(k) is a variational 
parameter. Magnons or spinon-spinon bound states are 
obtained from the holon-spinon bound states by creating 
an electron at the holon site, which forms a spin triplet 
with the spinon bound to it. 

Numerical comparisons of the trial wave functions for 
holes and magnons with the exact eigenstates are pre- 
sented in figure | and in tables 0-0. The hole trial 
wave functions for t±/J± = t/J = 1 (see table [n]) are 
excellent at stong coupling, and less accurate at isotropic 
coupling; for weak coupling, they are excellent only at 
momenta close to the one-hole ground states, as there 
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FIG. 4. Dispersions for a single (a) hole (holon-spinon 
bound state) and (b) magnon (spinon-spinon bound state) 
as predicted by the spin liquid proposed here in compar- 
ison with the exact dispersions for a 2 x 8 ladder with 
J± — J = t± — t — 1 (see tables |lj| and III) and periodic 
(or antiperiodic) boundary conditions. The dotted lines cor- 
respond to the individual + and — chirality trial wave func- 
tions (generated from the tight-binding configuration shown 
in figure ^b and its P or T conjugate, respectively) and the sta- 
tionary holetto-spinon bound state (according to figure ^a), 
as indicated. 



is a large amplitude to find a hole and a magnon rather 
than just a hole at other momenta. The holon-spinon 
bound state wave functions adequately describe the hole 
when t and J are comparable (see table IV); for t <C J, 
the holetto-spinon bound state is more appropriate, while 
holes with sufficiently large t are detrimental to antiferro- 
magnetic correlations and eventually destabilize the spin 
liquid pc| ] . The trial wave functions for the magnons are 
generally satisfactory (see table 111). 

6. The P and T violation of the localized holon-spinon 
bound states, or the appearance of a chirality quantum 
number, is a physical property of the system; any real 
trial wave function for the localized hole would yield a 
dispersion Ek oc cos(fc.j;), and thus be inconsistent with 
the dispersion obtained by exact diagonalization (bottom 
curve in figure ^a) . The chirality quantum number is a 
manifestation of the fractional statistics [^lj of spinons 
and holons jl8|,^2| in dimensions greater than one; it de- 
termines the sign of the statistical phases acquired as 
they encircle each other. 

I am deeply grateful to Bob Laughlin and Sudip 
Chakravarty for many illuminating discussions. This 
work was supported through NSF grant No. DMR-95- 
21888. Additional support was provided by the NSF 
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TABLE II. Energy expectation values including individual 
contributions to the kinetic energy from chains and rungs and 
overlaps for the trial wavefunctions describing holon-spinon 
bound states in comparison with the exact one hole eigen- 
states for a periodic 2x8 ladder with t±/J±=t/J — l for 
three ratios J±/J. The energies are in units of max(Jx,J); 



the transverse momentum is always k y 
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state corresponds to a hole (holon-spinon bound state) plus 
a magnon (spinon-spinon bound state). The trial wave func- 
tions are given by IV'triai) = Af(\ip+) + a|i/;_)), where is 
the trial wave function constructed with the flux configura- 
tion shown in figure ^b and \ip~) its P or T conjugate. The 
boundary phase for the flux band structure before projection 
has been for the chain containing the holon, and tt for the 
other chain. 
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TABLE III. As in table [n], but now for spinon-spinon 
bound states (magnons) with k y = n obtained by creation 
of an electron at the holon site. For k x — 0, spinetto-spinon 
bound states yield better trial wave functions (the energy is 
only 0.4 % off at J ± /J = 0.2, and 2.7 % off at J ± /J = 1). 



t/J 


Etot 


% 


over- 


E t 




E, 




M 




exact 


trial 


off 


lap 


exact 


trial 


exact 


trial 







-8.18 


-7.88 


3.7 


0.888 


0.00 


0.00 


0.00 


0.00 


1.00 






-8.00 


2.2 


0.938 




0.00 




0.00 


* 


0.2 


-8.37 


-8.14 


2.8 


0.927 


-0.17 


-0.14 


-0.10 


-0.16 


0.34 






-8.12 


2.9 


0.923 




0.00 




-0.13 


* 


0.5 


-8.87 


-8.63 


2.7 


0.929 


-0.65 


-0.53 


-0.28 


-0.33 


0.12 






-8.32 


6.3 


0.862 




0.00 




-0.32 


* 


1 


-9.89 


-9.52 


3.7 


0.894 


-1.54 


-1.23 


-0.59 


-0.57 


0.02 


2 


-12.1 


-11.3 


6.5 


0.797 


-3.36 


-2.63 


-1.31 


-1.02 


0.04 


5 


-19.5 


-16.8 


14 


0.585 


-8.86 


-6.84 


-3.76 


-2.33 


0.08 



TABLE IV. As in table but now with J±/J = t±/t = 1 
for different ratios t/J. The momentum is (k x , k y ) = (tt/2, 0), 
which corresponds to the ground state of the 2x8 ladder 
with periodic boundary conditions. For t/J < 0.5, data for 
holetto-spinon bound states are shown as well, marked with *; 
these adequately describe stationary holes (t — 0). 
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